Survival Probabilities of History-Dependent Random Walks 
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We analyze the dynamics of random walks with long-term memory (binary chains with long-range 
correlations) in the presence of an absorbing boundary. An analytically solvable model is presented, 
in which a dynamical phase-transition occurs when the correlation strength parameter y, reaches a 
critical value fi c . For strong positive correlations, fj, > fi c , the survival probability is asymptotically 
finite, whereas for fi < fj, c it decays as a power-law in time (chain length). 



Dynamical systems with long-range spatial or tempo- 
ral correlations are attracting considerable interest across 
many disciplines, lending applications to physical, biolog- 
ical, social, and economical sciences (see, e.g., [1-7] and 
references therein). Such systems are often analyzed us- 
ing random walks, a fundamental concept of statistical 
physics. Random walks lend applications to numerous 
scientific fields (see, e.g., [8-13] and references therein). 
In particular, random walks in the presence of absorbing 
traps (boundaries) have been studied in recent years as 
models for various systems such as absorbing-state phase 
transitions [14,15], polymer adsorption [16], granular seg- 
regation [17], the spreading of an epidemic [18], and in 
the context of complex adaptive systems [19,20]. In this 
work we analytically study random walks with an absorb- 
ing boundary, in which the jump probability is history- 
dependent, resulting in long-range correlations. 

The statistical properties of data such as DNA strings, 
written texts and financial data (e.g., stock market 
quotes), are known to significantly deviate from those of 
purely random sequences [21,22,7]. Such systems may be 
studied by mapping them onto a correlated sequence of 
symbols. Although the nature of the resulting sequence 
may depend upon the choice of mapping (see, e.g., [23]), 
the essential statistical properties of the original system 
are often preserved. By choosing a mapping of these sys- 
tems onto two symbols [5] , the problem is reduced to the 
exploration of correlated binary chains, which are equiv- 
alent to one dimensional random walks with a constant 
step size. These binary chains have long-range correla- 
tions and often exhibit a super-diffusive nature, in which 
the variance grows asymptotically faster than the string 
length. 

The preceding discussion motivates a study of random 
walks with a history-dependent jump probability, both 
with and without an absorbing boundary, as means of 
facilitating our understanding of systems with long-range 
correlation. In Ref. [7] we presented a model for a history- 
dependent random walk. Although simple, this model 
features a dynamical phase transition between normal 
diffusion and super-diffusion as the correlation strength 
parameter reaches a critical value. In this work we an- 
alyze random walks with long-range correlations in the 



presence of an absorbing boundary. 

We begin by introducing a simple model that incorpo- 
rates long-range correlations into an otherwise random 
sequence. Consider a discrete string of binary symbols, 
a,; 6 {0, 1}, in which the conditional probability of a 
given symbol (say, zero) occurring at the position L + 1 
is history-dependent, and is given by 



. 1/ L-2k 
p( k ,L) = -[l-n TTI - 



(1) 



where k is the number of such symbols (zeros) appear- 
ing in the preceding L bits. The correlation parameter 
/i, where — 1 < p, < 1, determines the strength of corre- 
lations in the system. The persistence condition // > 
implies that a given symbol in the preceding sequence 
promotes the birth of a new identical symbol. In the 
anti-persistence regime \i < 0, on the other hand, each 
symbol inhibits the appearance of a new identical sym- 
bol. The parameter Lq > 1 is a constant transient time. 
For L <C Lq the sequence is approximately random (un- 
corrected) , whereas for L ^> Lq the effect of correlations 
takes over [24]. 

In this model, the conditional probability p(k, L; /i, Lq) 
depends on the number of zeros (or unities) in the preced- 
ing bits, and is independent of their arrangement. This 
allows one to obtain an analytical description of the sys- 
tem's dynamical behavior. As demonstrated in [7], this 
two-parameter model provides a good description of the 
observed statistical properties of various systems such as 
coarse-grained DNA strings, written texts, and financial 
data, when mapped onto a binary chain. 

The probability P(k, L + 1) of finding k identical sym- 
bols (say, zeroes) in a sequence of length L + 1 follows 
from the evolution equation 



P(k, L + 1) = p(k - l,L)P(k - l,L) 
+[l-p(k,L)]P(k,L) . 



(2) 



Crossing to the continuous limit, one obtains the Fokker- 
Planck diffusion equation for the correlated process, 
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where x = 2k — L is the distance from the origin in the 
corresponding random walk, and we have neglected high 
order terms which are irrelevant for 1 <C Lq <C L. Solu- 
tions of Eq. (3) under the initial condition P(x, L = 0) = 
5(x) were given in [7]. 

We introduce an absorbing boundary at x = 0, by 
imposing the boundary condition P(x < 0, L > 0) = 
0. The evolution equation (3) along with this boundary 
condition and the initial condition P(x,L = 0) = S(x), 
has a solution in the form 



P (X, L) OC ^^3/2 x@ ( x ) ex P 



2D(L) 



(4) 



where 8 is the Heaviside step function, and D(L) is given 
by 



D(L;/j,,L ) 



L + Ln[ 
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L + Lr 



(5) 



Equation (5) breaks down in the special case ii 
where ii c = 1/2. In this case, D(L) is given by 



D(L\ n = fx c , L ) = (L + L ) In 



L + L 

Ln 



(6) 



The first two moments of the distribution function of 
the survived walkers are given by (x) 2 = ttD(L)/2 and 
(x 2 ) = 2D(L). The variance of the probability distribu- 
tion P(x, L) thus equals 



V(L;n,L ) = ^-^D(L) 



(7) 



This result implies that for fi < li c , the asymptotic vari- 
ance scales linearly with the string length, whereas for 
/i > /i c it scales as V oc L 2p . Hence, a history-dependent 
sequence with strong positive correlations (fx > ji c ) is 
characterized by a super-diffusion phase in which the 
variance grows asymptotically faster than L, both with- 
out [7] and with an absorbing boundary. 

The survival probability S(L) = f Q P(x,t)dx of the 
walkers is given, for fi fj, c , by 



L + Lr 



1-2(1 



-1/2 



whereas for /j = /i c we find 



S(L;[i c ) oc In 



L + L c 
Ln 



(8) 



(9) 



The survival probability thus changes its asymptotic 
(L 3> Lq) behavior at the phase transition value 
fx c = 1/2, and one identifies three qualitatively different 
regimes, 

( n< n c ; 

5(L»L )oc <^ l n -i(L/L ) fi = ^ ; (10) 
t const. [X > fi c . 



The normalization of the survival probability is sensi- 
tive to the discrete initial conditions. Since the chain is 
nearly random for L < L ! the normalization may be 
approximated by equating S(L), for 1 <C L -C Lq, to the 
survival probability of a purely random, continuous walk. 
We thus find 



S(l< L < L Q ) ~ erf 



(Vol/ 



(11) 



where xq is the distance between the absorbing boundary 
and the origin. 

In order to confirm the analytical results, we perform 
numerical simulations of (discrete) binary sequences. 
Figure 1 displays the resulting survival probability S(L) 
of correlated strings with various values of the correlation 
parameter fi. We find an excellent agreement between 
the analytical results [Eqs. (8), (9) and (11)] and the 
numerical ones. 




FIG. 1. Survival probability S(L) as a function of sequence 
length L. Shown are numerical (solid lines) and analytical 
(dashed lines) results for p = -0.8, -0.4, 0, 0.25, 0.5, and 0.8 
(from bottom to top). Here, xq — (half the walkers survive 
the first step) and Lq = 100. The analytical results agree 
with the numerical ones to better than 1%. 

The preceding results, starting from the distribution 
function in Eq. (4), are valid only when the absorbing 
boundary is placed at the origin. One would like to gen- 
eralize the results for boundaries located at arbitrary lo- 
cations xq, under the boundary condition P(x < xq, L > 
0) = 0. Note that the diffusion equation (3) is not in- 
variant under the translation x — * x + d, where d is a 
constant. The generalization of our solution for xq ^ 
is therefore non-trivial. Nevertheless, we have verified 
numerically that the asymptotic behavior of the survival 
probability given by Eq. (10) remains valid for arbitrary 
values of Xq, in all three regimes. The similar asymp- 
totic behavior for different choices of xq is demonstrated 
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numerically in Fig. 2. 




FIG. 2. Survival probability S as a function of string length 
L for various boundary locations xo- We present results for 
xo = 0, —2, —4, and —6 (from bottom to top), with fi — 0.2 
and L = 100. 

The above results can be readily generalized for a bi- 
ased random walk with a moving boundary. For example, 
for the biased jump probability 



p(k,L) 



1 + q — fi- 



L~2k 



(12) 



L + L 0/ 

the above results will hold if we apply the transformation 
x — > x — x c (L), with 



x c (L) = q 



L + L 



1 



/' 



L + L 



(13) 



In summary, in this work we have analyzed the dy- 
namics of random walks in which the jump probabili- 
ties are history-dependent, in the presence of an absorb- 
ing boundary. Using an analytically solvable model, we 
identify a dynamical phase transition characterizing the 
system's global behavior. For small values of the corre- 
lation strength (/i < fi c ) the survival probability decays 
as S cx L^^"^^, whereas for /i > fi c the system is char- 
acterized by finite asymptotic survival probabilities. 
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